Abstract 

This paper presents an algorithmic method to study structural proper- 
ties of nonlinear control systems in dependence of parameters. The result 
consists of a description of parameter configurations which cause differ- 
ent control-theoretic behaviour of the system (in terms of observability, 
flatness, etc.). The constructive symbolic method is based on the differ- 
ential Thomas decomposition into disjoint simple systems, in particular 
its elimination properties. 

1 Introduction 



Symbolic computation allows to study many structural aspects of control sys- 
tems, e.g., controllability, observability, input-output behaviour, etc. In contrast 
to a numerical treatment, the dependence of the results on parameters occurring 
in the system is accessible to symbolic methods. 

An algebraic approach for treating no nlinear control s y stems has been de 



velop e d during the last decades, cf., e.g., iFliess and GladI (|1993[ ). iDiod (|1991 



1992f ). Pommaret ( 2001 ). and the references therein. In particular, the notion of 
flatness has been studied extensively and ha s been applied to many interesting 



control problems (cf., e.g., 
the characteristic set method (cf. iKolchi: 
Grobner algorithm (cf. iBoulier et al. ( 200] 



Fliess et al] ( 19951) ). The approach of Diop builds on 



Wul (I2OOOI) ). The Rosenfeld- 



can be used to perform the relevant 



computations effectively; cf. also Wand ( 2001 ) for alternative approaches. 

Dep endence of control systems on parameters has bee n examined, in partic- 



ular, in lPommaret and Quadrat) (jl997l ). iPommaretl (|200ll) . In the case of linear 



systems, stratific ations of the space of parameter values have been studied using 
Grobner bases in iLevandovskvv and Zera (j2007l) . 

In the 1930s the American mathematician J. M. Thomas designed an al- 
gorithm which decomposes a polynomially nonlinear system of partial differen- 
tial equations into so-called simple systems. The algorithm uses, in contrast 
to the characteristic set me thod, ine q uation s to provide a disjoint decom posi- 
tion of the soluti on set (cf. iThomas (Il937 )'l. It precedes Kolchin ( 1973 ) and 
Seidenbere ( 19561 ) (building on Ritt (jl95C )). Recently a new algorithn iic ap- 
proach to the Thomas decomposition meth od has been developed (cf. iGerdtl 
(|2008|) : iBachler et all (l2012h: lRobert3 (|2012l )). also using ideas of M. Janet (cf. 
Janet! (Il929l) ): cf. also IWana (|1998l ) for an earlier implementation of the alge- 
braic part. 

So far the dependence of nonlinear control systems on parameters has not 
been studied by such a rigorous method as the Thomas decomposition. This 
paper demonstrates how the Thomas decomposition method can be applied 
in this context. In particular, the Thomas decomposition can detect certain 
structural properties of control systems by performing elimination and it can 
separate singular cases of behaviour in control systems from the generic case 
due to splitting into disjoint solution sets. 

In Section [2] we sketch the Thomas decomposition method for algebraic and 
differential systems. The algorithm for the differential case builds on the al- 
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gebraic part. Section [3] describes how the differential Thomas decomposition 
can be used to solve elimination problems that occur in our study of nonlinear 
control systems. In Section|3]we recall the notions of observability and flatness, 
which are addressed in the examples in Section [S] using our Maple implementa- 
tion. 

2 Thomas Decomposition 

2.1 Simple Algebraic Systems 

Let i^T be a field of characteristic zero. We denote by R the polynomial algebra 
K[xi, . . . , Xn] and we fix a total ordering < on {xi, . . . , Xn}- Then the greatest 
variable with respect to < which occurs in a non-constant polynomial p £ R is 
called the leader oi p and is denoted by ld(p). 

In what follows, we regard every p d R \ K as a, polynomial in ld{p) with 
coefficients in K[xi \ 1 < i < n, Xi < ld(p)], and, recursively, each coefficient 
as a polynomial in its leader, etc. The coefficient of the highest power of \d{p) 
in p is called the initial of p, denoted by init(p). Moreover, we denote the 
discriminant oi p (as a polynomial in ld(p)) by disc(p). Both init(p) and disc(p) 
are elements of the polynomial ring K[xi \ 1 < i < n, Xi < ld{p)]. 

Let K be an algebraic closure of K. In this subsection we consider (finite) 
algebraic systems (of equations and inequations, defined over K) 

5 - {pi - 0, . . . , = 0, <zi ^ 0, . . . , 7^ 0}, (A) 

where pi, qj £ R, s, t E Z>o. The set of solutions of S (in K ) is defined as 

Sol^(5) = {a e i?" I p,{a) = 0, qj{a) ^ for aU i,j}. 

Finally, for < fc < n, let tt^ : k" — > k" ^ be the projection onto the 
coordinate subspace of dimension n — k whose coordinates are ranked lowest 
with respect to <, i.e., 7rfc(ai, . . . , a„) = (oi^. | 1 < j < n — k), where 1 < «i < 
. . . < in~k < and Xi^, . . . , Xi^_^. are the smallest n—k elements of {xi, . . . , Xn} 
with respect to <. 

Definition 1. An algebraic system S as in is simple (with respect to <) if 
the following conditions hold. 

1. Pi ^ K, Qj ^ K for all i and j. 

2. |{ld(pi),...,ld(p,),ld(gi),...,ld(qO}| + 

3. For every r S {pi, . . . qi, . . . , }, if ld(r) is the fc-th greatest vari- 
able w.r.t. <, then neither init(r) = nor disc(r) — has a solution in 
^fe(Sol^(5)). 
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Remark 2. Condition (2) implies that the leaders of the left hand sides of the 
equations and inequations in a sim ple algebraic systeni S are pai r wise d i stinct. 



S is a triangular set (cf., e .g., lAubrv et all (|l999l) : iHubertI (|2003l ): IWand 



( 2001 ) ; cf . Lemaire et al. ( 20051) for a related implementation) . 

The meaning of condition (3) is that for every 1 < A: < n and each (ai, . . . , Qn-k) 
7rfe(Soly(5)) there exists aGK such that (oi, . . . , a„_/c, a) e 7rfe_i(Sol-^(S')), 
and the number of possible values for a is either finite or co-finite, determined 
by the degre e of the leader in the corresponding equation or inequation in 5*, if 
any (cf. also IPleskenI |2009) ). 



Definition 3. Let S be an algebraic system as in A Thomas decomposition 
of S (with respect to <) consists of finitely many simple algebraic systems 
5*1, . . . ,Sk such that Sol-j^ (5*) is the disjoint union of Sol-f^(5'i), . . . , Solj^^Sk)- 

Remark 4. Euclidean pseudo-division and splitting of systems allow to com- 
pute a Thomas decomposition for any given (finite) algebraic system (defined 
over a computable field K of characteristic zero) in finitely many steps. 

In each round the algorithm chooses a system S from a set of systems to be 
processed. If S contains an equation whose left hand side is a non-zero constant 
or an inequation with zero left hand side, then S is discarded because its solution 
set is empty. Moreover, we assume that equations = and inequations k ^ 0, 
where k G K \ {0} are removed from each system, so that now the leader of 
every equation and inequation is well-defined. 

The algorithm applies Euclidean pseudo-division to each pair pi — 0, p2 — 
of distinct equations in S with Id(pi) — ld(p2) =: x] i.e., if deg^(pi) > deg^(p2)) 
then polynomial division is performed on c-pi and p2, where c is a suitable power 
of init(p2)- In order not to change the solution set, when = is replaced 
with the result of the polynomial division, non- vanishing of init(p2) on the set 
of solutions of S is assumed. To this end, a preparatory step splits a system 
into two and adds init(p2) 7^ and init(p2) — 0, respectively, if necessary. 

Each pair qi ^ 0, q2 ^ of distinct inequations in S with \d{qi) = ld(g2) 
is replaced with lcm(gi,g2) 7^ 0. The computation of the least common multi- 
ple depends on case splittings according to vanishing of initials, similar to the 
previous remarks. 

For each pair p — 0,qj^OmS with \d{p) ~ ld((7), the algorithm computes 
r := gcd(p, q) using Euclidean pseudo-division. This process again relies on 
non-vanishing initials, which is ensured by distinguishing cases appropriately 
in advance. If p divides q, then S has no solution and is discarded. If r is a 
non-zero constant, then (7 7^ is removed from S. Else, p — and g 7^ are 
replaced with p/r = and q/r 7^ 0, respectively. 

Finally, there is some flexibility when to take care of non-vanishing discrim- 
inants. Since K is of characteristic zero, this property is established for any 
equation r = or inequation r 7^ by examining g := gcd(r, ^jg^), again with 
appropriate case distinctions, and thus determining the square-free part r/g. 
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We refer to iBachler et alj (I2OI2I): IPleskeij (l2009h: iRobertj (|2012l) for more 
information and to Bachler and Lange-HegermannI ( 201 2! ) for an implementation 
in Maple. In practice, we apply subresultants for the computation of gcds and 
the related case distinctions. 



2.2 Simple Differential Systems 

Let if be a differential field of characteristic zero with pairwise commuting 
derivations 61, . . . , (5„, i.e., A" is a field of characteristic zero, and each Si is 
a map K ~> K satisfying Si{ki + — Si{ki) + 6i{k2) and the Leibniz rule 
6i{ki k2) = ki 5i{k2) + Si{ki) ^2 for aU fci, ^2 € K. 

In this subsection we denote by R the differential polynomial ring K{ui, . . . , u^} 
in the differential indeterminates ui, . . . , u„i with pairwise commuting deriva- 
tions di, dn whose restrictions to K are Si, (5„, i.e., R = K[{ui)j \ 
1 < i < TO, J G (Z>o)"] is the polynomial ring in the algebraically indepen- 
dent variables {ui)j, and dj : R R is the derivation defined by extending 
dj{{ui)j) := {ui)j+ij additively to R such that it satisfies the Leibniz rule on 
R and restricts to Sj on K. Here Ij is the multi-index of length n whose i-th 
entry is 1 if i = j and otherwise. 

Each variable {ui)j is thought of as representing the partial derivative, cor- 
responding to the multi-index J, of a smooth (or rather analytic) /^-valued 
function of n arguments. When dealing with differential polynomials p € R 
algorithmically, we compare terms with respect to a ranking, which is defined 
next. We agree that applying derivations should make terms larger with respect 
to the ranking. This is taken into account as follows. 

Definition 5. A ranking on i? is a total ordering < on {{ui)j | 1 < i < to,. J € 
(Z>o)"} such that Ui < djUi for all i and j, and such that (uijji < (wi2)j2 
implies dj{ui^)j^ < dj{ui^)j^ for ah j, ii, 12, Ji, J2- 

Every ranking is a well-ordering, i.e., there exist no infinitely decreasing 
sequences of variables {ui)j. 

We define := d-^K..d^^, J e (Z>o)", and write | J| Ji + . . . + J„ for 
the length of the multi-index J. 

Example 6. The degree-reverse lexicographical ranking on the differential poly- 
nomial ring K{u} (i.e., to = 1) with pairwise commuting derivations di, . . . , 9„ 
is defined by 

r |J| < |J'| or (l^l = k'l and J 7^ J' and 
uj < uj, :^ I ^ J, i = max{l < k < n \ Jk ^ J^}). 

In this example, we have 9„m < 9„_im < . . . < diu. There are in fact nl different 
degree-reverse lexicographical rankings according to the ordering of the diU. 

In this subsection we consider (finite) differential systems (i.e. systems of 
differential equations and inequations) 

S = {pi = 0, ...,p, = 0,qi^O, ...,qt^O}, (D) 
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where pi, qj € R, s, t € Z>o. 

In what foUows, we assume that a ranking < on R is fixed. Carrying over 
the concepts of Subsection 12. 1[ the leader \d{p) of p G R\K is defined to be the 
greatest variable with respect to < which occurs in p. The initial of p, denoted 
by init(p), is the coefBcient of the highest power of ld(p) in p, and the separant 
of p, denoted by sep(p), is the formal partial derivative of p with respect to 
ld{p). The initial of p is an element of the polynomial ring K[{ui)j | 1 < i < 
m, Je (Z>o)", K)j <ld(p)]- 

In order to ensure formal integrability for the kind of differ ential s ystem s 
we are heading for we a pply the concept of Janet division (cf. IJanetj ( 19291 ): 



Gerdt and Blinkovl (|l998f )). which restricts the usual divisibility relation on the 



free commutative monoid O generated by 9i , . . . , 9„ . 

Definition 7. Let M C (Z>o)" be finite. Janet division associates with each 
/ G A/ a partition of into the subsets of admissible and non- 

admissible derivations as follows. Let / = G M. Then di is ad- 

missible for / if and only if 

h ^ max{/; !(/(,..., O G A/, I'j = Ij for all 1 < j < i}; 

otherwise di is non-admissible for /. We denote by M) the set of derivations 
that are admissible for / (with respect to M), Ji{I, M) := {di, . . . , 9„}\/i(/, M), 
and 

e(/, M) ■= {d-^ I J G (Z>o)", Jj = if dj G 7I(/, M)}. 

Janet division defines for each with / G A/ a cone of multiples Q{I, M)d^ 
such that each two cones are disjoint. By enlarging M appropriately, these 
cones cover the set of all multiples of 9^, / G M . 

Definition 8. A finite subset M of (Z>o)" is said to be Janet- complete if 



Definition 9. Let pi, . . . , pr G R\K , ld{pj) = d^^Ui. , be such that each Mk = 
{Ij I 1 !i J < *j = k} is Janet-complete. A differential polynomial p € R 
is Janet-reduced modulo pi, . . . , pr if, for all 1 < j < r, deg^(p) < deg^(pj), 
where x := Id(pj), and no variable d''\d{pj), where d'' G Q{Ij,Mi.), \ J\ > 0, 
occurs in p. 



Remark 10. Any p E R can be transformed into a differential polynomial 
which is Janet- reduced modulo pi, . . . , p^hy applying Euclidean pseudo-division 
modulo pi, . . . ,pr and their derivatives repeatedly. The obvious strategy is to 
first eliminate the greatest variable d'^\d{pj) with respect to < that occurs in p 
and to proceed to lower variables. If d'^\d{pj) is a proper derivative of \d{pj), 
then, before reducing modulo pj, the pseudo-division multiplies p with sep(pj) 
(which, in fact, is the coefficient of the leader of any proper derivative of pj). If 
I J| = 0, then the pseudo-division multiplies p with init(pj). As in Subsection l2.1l 
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cases of vanishing initials or separants have to be examined separately, which 
results into splittings of differential systems. Under the assumption of non- 
vanishing initials and separants we denote the result of the pseudo-reduction by 

Definition 11. In the situation of Definition^ {pi, . . . ,pr} is passive if, for all 
1 < j < r, NF{dpj,{pi, . . . = for aU d e Jl{Ij,M,^). 

Definition 12. We call a differential system S as in ([D|) simple (with respect 
to <) if the following conditions hold. 

1. 5 is simple as an algebraic system (in the finitely many variables {ui)j 
which occur in S, totally ordered by the ranking <). 

2. {pi, . . . ,ps} is passive. 

3. qi, . . . , qt are Janet-reduced modulo pi, . . . , p^. 

From now on we assume that K is a differential field of complex meromor- 
phic functions on a connected open subset D, of C" with coordinates zi, . . . , 
Zn, and that the derivations di, (5„ are defined by partial differentiation 
with respect to zi, . . . , z„, respectively. Differential equations p — (and their 
derivatives) translate into algebraic equations for the Taylor coefficients of a 
power series expansion of ui, . . . , Um around an arbitrary point in fl. A dif- 
ferential inequation q ^ can be interpreted as the disjunction of algebraic 
inequations for all Taylor coefficients of the analytic expression that is obtained 
from q by substitution of power series expansions for ui, . . . , Um- 

In considering complex analytic functions on fl as solutions to differential 
systems, we assume that is chosen appropriately with regard to given dif- 
ferential systems, more precisely, that for every problem instance treated by 
the Thomas decomposition method, for each of the resulting simple differential 
systems there exists an analytic solution on ft. Which subsets il of C" are 
appropriate can often be decided only after the formal treatment of the given 
differential systems by the Thomas algorithm. Questions concerning the radius 
of convergence of analytic solutions are ignored here. 

For any differential system S as in ([D|) we denote by Solsi(S') the set of 
complex analytic functions / : J7 -^• C satisfying Pi{f) — and qj{f) ^ for all 
i and j. 

Definition 13. Let S* be a differential system as in (|D|. We call a family of 
finitely many simple differential systems Si, . . . , Sk such that 8010(5) is the dis- 
joint union of the solution sets 8010(5*1), . . . , 8olo(5fc) a Thomas decomposition 
of S (with respect to <). 
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Remark 14. Given any (finite) differential system (with coefficients in a com- 
putable differential subfield of i^), a Thomas decomposition into simple dif- 
ferential systems can be computed in finitely many steps by a process which 
interweaves the algebraic Thomas algorithm (cf. Remark 2]) and Janet reduc- 
tion (cf. Remark [TU|) . 



Theorem 15. fcf. \Roberd tOlj )) Let S as in W\) be a simple differential 
system, E the differential ideal of R generated by pi, . . . ,ps, and let q be the 
product of all mit{pi) , sep(pi). Then 

E : := {p e R \ q'' ■ p e E for some r G Z>q} 

is a radical differential ideal, which consists of all differential polynomials in 
R vanishing on 8010(5*). Given p G R, we have p G E : q°° if and only if 
NF(p, {pi,...,p,}) = 0. 

We refer to iBachler et all (|2012|): iGerdtl (l2008l): iRobertj (|2012l ) for more 
information and to Bachler and Lange-Hegermann (|2012l ) for an implementation 
in Maple. 



3 Elimination 

We continue to use the same notation as in the previous section. Our objective 
is to perform a projection of the solution set of a differential system onto the 
space which is addressed by only certain of the components of the solution 
tuples. In other words, we would like to determine all differential consequences 
of the given system involving selected differential indeterminates only. 

Definition 16. Let Bi, . . . , Bk form a partition of the set of differential in- 
determinates {ui, . . . , Ufn}- The block ranking on R with blocks Bi, . . . , Bk is 
defined as follows, where G Bj^, G Bj^, Ji, J2 G (Z>o)": 

d-'^u,^ < d-'^u,^ :^ ii > js or (ji = and d'^^ < d'^^). 

The comparison of d^^ and d^'^ is defined by a choice of a degree-reverse 
lexicographical ordering as in Example [51 Such a ranking is said to satisfy 
> S2 > . . . > B/c- 

Using the above notation, for any 1 < i < /c, we define K{Bi, . . . , Bk] 
K{u \uG BiU...UBk}CR. 



Proposition 17. fcf. iRobert^ \201^ )) In the situation of Theorem ] 151 suppose 



that < is a block ranking with blocks Bi, . . . , Bk- For 1 < i < fc, let Ei be the 
differential ideal of K{Bi, . . . , B^} generated by {pi, . . . ,ps} H K{Bi, . . . , Bk}, 
and let qi be the product of all initials and separants of the elements in this 
intersection. Then, for every 1 < i < k, 

iE:q^)nK{B„...,Bk}^Er.qr- 
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Hence, computing a Thomas decomposition with respect to a block ranking 
enables us to extract generating sets for the differential consequences satisfied 
by the projection of the solution set of the given differential system. 

For any differential system S let (resp. S^) denote the set of left hand 
sides of equations (resp. inequations) in S. 



Corollary 18. fcf. \Roberd kOlj )) Let S he a (finite, not necessarily simple) 



differential system, and let Si, . . . , Sr be a Thomas decomposition of S with 
respect to a block ranking with blocks Bi, . . . , ■ Let E be the differential 
ideal of R generated by and define the product q of all elements of . Let 
i e {1, . . . , fc}. For 1 < j < r, let Ej be the differential ideal of K{Bi, . . . , Bk} 
generated by SJ" D K{Bi, . . . , B^}, and let qj be the product of all initials and 
separants of the elements in this intersection. Then, 

^/ET^nK{B„...,Bk}^{Ei:q^)n...niEr:qn. 



4 Systems Theory 

In this section we adapt well-known co ncept s of no n hnear control theory t o 
our fram ework (cf. e.g., iGlad' (l9M); IPiopI (Il992h : iFliess and GladI 



Fliess et a l. (19951: iPomm arct (2001D and the references therein). 

Let R be the differential polynomial ring K{ui, . . . ,Um} in the differential 
indeterminates ui, . . . , Um with pairwise commuting derivations di, . . . , (?„ as 
in Subsection 12.21 

We assume that a (nonlinear) control system is given by a simple differential 
system S as in ([D]). Let E be the differential ideal of R generated hy pi, ... ,ps, 
and let q be the product of all init(pi), sep(pi). Let U := {ui, . . . ,Um}- (No 
distinction is made a priori between input, output, state variables, etc.) 

Definition 19. Let x G [/ and Y (- U \ {x}. Then x is observable with respect 
to Y if there exists p G {E : q°°) \ {0} such that p is a polynomial in x (not 
involving any proper derivative of x) with coefficients in K{Y} and such that 
neither its leading coefficient nor || is in i? : q°°. 



Remark 20. Given the components of a solution to S corresponding to the 
variables in Y , the implicit function theorem yields the component correspond- 
ing to a;, if there exists a polynomial p as in Definition [121 If S is simple with 
respect to a block ranking satisfying U \ {Y U {x}) ^ {x} ^ F by Corol- 
lary [TH] there exists such a p G {E : q°°) \ {0} if and only if such a p exists in 
n K{Y,x}. Otherwise, one has to compute a Thomas decomposition of S 
with respect to such a ranking and inspect each simple system. 



Definition 21. Let F C t/. Then Y is called a/laf output if {E : q°")r\K{Y} = 
{0} and every x E U \ Y is observable with respect to Y. The control system 
given by S is said to be flat if a flat output exists. 
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Remark 22. Whereas, again, computation of a differential Thomas decompo- 
sition of S with respect to a block ranking satisfying U\Y ^ Y allows to decide 
whether y is a flat output, deciding whether a given nonlinear control system 
is flat is a difficult problem in general. 

Many further applications of the Thomas decomposition method to systems 
theory (e.g., computation of the input-output behaviour from a state space 
representation, parameter identification, realization, inversion) can be realized 
and will be studied in the future (cf. also Diop ( 1992I )). 



5 Applications 

In this section we demonstrate the Thomas decomposition method on two ex- 
amples. 

Example 23. As an ap plication we consider the mode l of a continuous stirred- 



tank reactor taken from iKwakernaak and SivanI (1972) 



This model describes a tank with a dissolved material of concentration c, 
which is assumed to be the same everywhere in the tank due to stirring. Two 
input feeds with flow rates Fi and F2 feed this material into the tank with 
constant concentrations Ci and C2, respectively. There exists an outward flow 
with a flow rate proportional to the square root of the volume V of liquid in the 
tank. The system is modelled by the two differential equations 

V{t)^Fiit)+F2{t)-ky/V(tj 



C{t) V{t) = Ci Fi{t) + C2 F2{t) - c{t) k v/^O , 

for an experimental constant k. 

In the following we want to describe the properties of the system in depen- 
dence of the constants ci and C2 . The Thomas algorithm performs this analysis 
if we model these contants as functions ci{t) and C2(t) satisfying the differen- 
tial equations ci{t) — and C2{t) — 0. Additionally, the square root y/V{t) of 
V{t) appears in the equations. Our formalism cannot handle roots of functions 
directly; instead, we introduce \/V{t) as new differential indeterminate and 

substitute V{t) by (^\/V{t)^ . To exclude trivial cases, we assume ci{t) 7^ 0, 

C2it) y^O, andVit) j^O. 

We apply our implementation (cf. iBachler and Lange-Hegermann ( 2012 )). 



The command ComputeRanking sets the ranking. Its first parameter is the 
list of independent variables and the second parameter is the list of differential 
indeterminates; a list of lists of differential indeterminates indicates a block 
ranking. We input a derivative of a differential indeterminate as the name of 
the differential indeterminate indexed by its order; for example sV[l] stands for 
■^y/V{t). The main command Dif f erentialThomasDecomposition computes 
a differential Thomas decomposition for a given list of equations and a list of 
inequations with respect to the fixed ranking. 
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In the following we compute a Thomas decomposition of the system using a 
ranking with {Fi, i^z} > {W, c} » {01,02}. 

> ivar:=[t] :dvar:=[[Fl,F2] , [sV,c] , [cl,c2]] : 

> ComputeRankingCivar ,dvar) ; 

> L: = [ 2*sV[l]*sV[0]-Fl[0]-F2[0]+k*sV[0] , 

> c [1] *sV [0] -2-c2 [0] *F2 [0] +c [0] *k*sV [0] 

> -cl [0] *F1 [0] +2*c [0] *sV [1] *sV [0] , 

> cl[l] , c2[l]] : 

> res : =Dif f erentialThomasDecomposition( 

> L, [sV[0] ,cl[0] ,c2[0]]) ; 

res := [DifferentialSy stem, DijferentialSy stem, 
DijferentialSystem] 

This yields a decomposition consisting of three simple systems. We print the 
first system; for better legibility we have underlined the leaders of equations. 

> subs(sV(t)=sqrt(V(t)) , 

> PrettyPrintDif f erentialSystem(res [1] ) ) ; 



[(c2(t)-cl(i))Fl(t)+ (^|c(t)) (v^)' 

+ m o2{t)) (^2j^VW) + fc) Vm = 0, 
(clW-c2(i))F2(t)+ (^|c(t)) {Vrnf 

+ m - oiit)) (^2 1 + fc) = 0, 



(1) 



(2) 



|cl(t).0,|c2(t).0, (3) 



V^^O, cl(t)-c2(t)^0 
o2{t) ^ 0, cl(<) ^ 0] 



The equations (P) and @ allow us to solve for Fi (t) and F2 (t) given any o{t) 
and V{t). Thus, we consider o{t) and V{t) a flat output of the system under the 
additional condition ci — C2 7^ on the constants. Note that the two equations 
ci{t) — and 62 (t) = in ([3]) just model the parameters ci and C2 as constants. 

The other two systems of this decomposition have the additional condition 
Cl = C2 on the constants. This condition prohibits to control the concentration 
in the tank as both input feeds are equivalent. Thus, these systems do not admit 
c{t) and V{t) as a flat output. 

Now we turn our attention to the observability of \/V{t) using Remark [20l 
Therefore, we choose a ranking with {a/V^} ^ {c, Fi, F2} ^ {ci, 02} ■ A Thomas 
decomposition with this ranking consists of seven systems: 
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> ivar:=[t] :dvar:=[[sV] , [c,Fl,F2] , [cl,c2]] : 

> ComputeRankingCivar ,dvar) ; 

> res : =Dif f erentialThomasDecomposition( 

> L, [sV[0] ,cl[0] ,c2[0]]) ; 

res :— [DifferentialSystem, DifferentialSystem, 
DifferentialSystem, DifferentialSystem 
DifferentialSystem, DifferentialSystem, 
DifferentialSystem] 

In the first two systems an equation with \/V{t) as leader appears and thus 



y/V{t) is observable. For the first system the condition on the parameters for 
observabihty is 

{c{t) - Ci)Fi{t) + {c{t) - C2)F2{t) ^ 0. 

The second system is not physically feasible as it involves negative input feeds 
due to i^2(0 and = -F2(t). 

The other five systems include an equation with ji\/vjt) as leader and 
thus \/V{t) is not observable. This follows for the first of these systems as it 
contains the Wronskian Fi{t)F2{t) — F2{t)Fi{t) = which makes the inputs 
linearly dependent. In the second system all three concentrations c{t), ci, and 
C2 are equal and constant. In the third system one input feed is zero and the 
concentration in the tank is equal to the concentration in the other input feed. 
The last two of these systems are not physically feasible because of negative 
values, as above. 



Example 24. Now we conside r a system of partial differential equations from 
IPommaret and Quadrat (1997). The system consists of three linear pdes in 



= £,i{xi,X2,X3) for i= 1,2,3 with a parametric function a{x2)- 

+ la{x2){V-asi)), (4) 

= ^a(x2)Ae.fe) + Ae3fe), (5) 

= -«(x2)^6fe) + - ^ (V • . (6) 

We model the parameter 0(2:2) as function a(a;i, a;2, 2:3) in three independent 
variables satisfying the differential equations 

d d 
-^a{xi,X2,X3) = and ■^-—a{xi,X2,X3) = 0. (7) 

Note that names for derivatives of differential indeterminates now involve multi- 
indices; for example xil [1,0,0] stands for gf^Cife)- We apply our implemen- 
tation with a block ranking {^1,^2,^3} ^ {a}- 
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> L:=[ -a[0,0,0]*xil[l,0,0]+xi3[l,0,0] 

> -a[0,l,0]*xi2[0,0,0]+(l/2)*a[0,0,0]* 

> (xil[l,0,0]+xi2[0,l,0]+xi3[0,0,l]) , 

> -a[0,0,0]*xil[0,l,0]+xi3[0,l,0] , 

> -a [0 , , 0] *xil [0 , , 1] + ( 1/2) *xi3 [0,0,1] 

> -(l/2)*xil[l,0,0]-(l/2)*xi2[0,l,0] , 

> a[l,0,0] ,a[0,0,l]] : 

> ivar : = [xl ,x2 ,x3] : dvar : = [ [xil ,xi2 ,xi3] , [a] ] : 

> ComputeRanking (ivar, dvar) ; 

> res:=DifferentialThomasDecomposition(L, [] ) ; 

res :— [DifferentialSy stem, DijferentialSy stem, 
DijferentialSystem] 

The Thomas decomposition yields three systems. The first system con- 
tains no additional condition for a{x2) except ([7]). In this generic case 
(O, and (O are not formally integrable leading to the compatibility condi- 
tion S,2{xi,X2,X3) = 0. The second system includes the additional condition 
^a(x2) = 0, which was calculated in IPommaret and Quadrat ( 1997 ) as the 



condition to ensure formally integrability of (|4]), ([5]), and (|6]). The third system 
is a special case of the second system with the new condition a{x2) = 0. 
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